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1 Usefulness of calculating Fourier series in the SEM 

In this note is presented a method, given nodal values on multidimensional 
nonconforming spectral elements, for calculating global Fourier- series coeffi- 
cients. This method is "exact" in that given the approximation inherent in 
the spectral-element method (SEM), no further approximation is introduced 
that exceeds computer round-off error. The method is very useful when the 
SEM has yielded an adaptive-mesh representation of a spatial function whose 
global Fourier spectrum must be examined, e.g., in dynamically adaptive fluid- 
dynamics simulations such as [Z|. 
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2 Derivation of an exact transform 



Suppose we have some functional problem in a spatial domain D := [— 7r,7r] rf 
(possibly including toroidal geometry) and use coordinate transformations 

d k from f G E := [-1, if to x G E k (1) 

to partition D = UjIILi E& by if elements : = $fc(Eo) with disjoint 1 interiors. 
Typically the SEM approximates the exact solution by its piecewise polyno- 
mial representation of degree P: 
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u cx (x ) « = ]T X] u T , k (l> T , k (x), (2) 



where J := {0, . . . P} d indexes the values u^ k : = u(x^ k ) and nodes sc^*. : = ^A:(^j*) 
mapped from the (i-dimensional Gauss-Lobatto-Legendre (GLL) quadrature 
nodes G [-1,1], 



(3) 



(j)yO'd k 1 {x) : X G E fc 
0, fG'Efc 
is the xp-interpolating piecewise-polynomial, 

M0--=IiMn and ^(O^EfepMO (4) 

a=l p=0 

are £y- and ^-interpolating polynomials, = WjL p (£j)/X]jLo w j'^>(£j') 2 * s 



a Legendre coefficient [e.g.,|j, (B.3.15)], yp + |L P (£) is the orthonormal Leg- 
endre polynomial of degree p on [—1,1] and Wj is the GLL quadrature weight. 
In many cases a physically interesting quantity is the global Fourier-series co- 
efficient u$ at integer wavenumber components q a , usually approximated by 



1 E k p| E fc / = if jfe ^ k' 
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M d -point trigonometric <i-cubature in such manner as 



1 



where 



Uq:-- 



K 



(o \d i n u ^ e iq ' x dv ^ = 



M d 



035* (^»») e -£tfl>f,ki 



(5) 
(6) 



dv(x) := Yia=i dx a is the volume differential and M := {1, . . . M} d indexes 
trigonometric nodes := (2m a /M — l)ir. Note whenever D is adaptively 
repartitioned there is an additional computation cost of 0(M d ) per node to 
use (2) to provide in (6) the values as well as a <i-cubature error 



[generalizing 
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theorem 4.7] 



fez d \{6} 

that in general converges no faster than 0(M~ 2 ), because C 1 discontinuities 
of (2) across element boundaries cause \uq-\ to decay only as 0{\q\~ 2 ). We 
discover a more accurate method by substituting (3) into (5) to yield 



-t&foo^Hx) dv(x) 



(i) 1 



(2n) d Jeo 

(4) 1 
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(2vr) rf 7e 



d p 

i a=l p=0 



) 









In many applications, especially when ^-structure rather than domain geom- 
etry is guiding the mesh adaption, each is a ^-parallelepiped with center 
c4 and d legs 2h^, so we have an affinity := + h^ • £, where fajj make 

up the columns of Then we obtain 



(2vr) c 



~^ fc n e 

a=l p=0 
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Finally, recalling the classical identity [e.g.,|l|, exercise 12.4.9] for the spherical 
Bessel function B p (r) of the first kind, 



B p (r) = ^£e--%(Od£, (7) 



we obtain ^ K<} = — 



Note that most expressions in (8) can be precomputed; objects that may vary 
during a dynamically adaptive computation, such as a k or h^, typically take 
values from a sparse set, e.g., a collection of powers of 2. The computation of 
(5) now incurs no additional error beyond that of (2). Also note, to generalize 
to the case P = P[? is straightforward. 



3 Accuracy of transform for ID & 2D test cases 

Equation (8) was implemented in MatLab® and tested using known results 
for (5). The most immediate test follows from (7), namely if* = L p (-/ir) = 
i _p Bp(7rg). In this case (5) was found to reproduce (7) to 12-16 digits for 
K — 1, P < 18, implying similar performance for any polynomial u(x) in this 
range. The next test was to put w ex (x) = sina;, or if* = (S q> i — <5 9j _i) /2i. Since 
this is not a polynomial we should expect at best to see algebraic convergence 
w.r.t. K in a uniform meshing a k = (k — l)h k — 7T,hk = 27r/K and exponential 
convergence w.r.t. P, as verified in Fig. 1. Note there is no need to test M ex (x) = 
sin rx for r > 1 because of scaling. 



We conclude by examining three 2D tests with adaptive meshing in the fashion 

n n 

of 5], using MatLab®. Fig. 2 confirms (5) in the case jy, (19)] 



u ex (x) = J2 e bl l <?1 l +b2 l g2 l +i ^, (9) 

72 



4 




6 



8 



14 



log 2 (element count K) 10 16 



Fig. 1. Surface plot (blue low to red high) of logio relative r.m.s. error in (5) for 
u ex (x) = sinx, vs log 2 K and P. 

where b a = — | and I = (pfi) = (-21) is a biperiodicity-preserving "rota- 
tion". As expected, the red curve (connecting the \u^\ peaks) shows a power- 
law decay in g^-space. Note, in this plot and those below the /-operation helps 
instigate mesh adaption but has the consequence of leaving q undersampled 
in Z 2 . In Fig. 3 is shown an initial condition P, (22)] 

{T x (0,f) := -tsinf-x (10) 

for the 2D Burgers eq. As expected, almost vanishes for (f 7^ ±Z. Finally, 



2, (2.5)] to 2D is 



at time t = 1.6037/ 7r|Z| the analytic solution generalizing 
shown in Fig. 4. As expected for the nearly C°-discontinuous fronts _L I seen 
at left, decays slightly faster than (9(|g| _1 ) but only for wavevectors q\\l 



(red curve). 
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Fig. 2. Left, u (9) over the spatial x domain, increasing from blue to red; yellow 
lines indicate element boundaries, black lines show nodes xj^ with P = 5. Right, 
surface plot of {uA from (5) vs q 1 and q 2 . 




Fig. 3. As in Fig. 2 but for the t = state given by (10), in K = 2 6 elements. 
[2] C. Basdevant, M. Deville, P. Haldenwang, J.M. Lacroix, J. Ouazzani, R. 

Peyret & P. Orlandi, 1986: Spectral and finite difference solutions of the 

Burgers equation. Computers & Fluids, 14, 23-41. 

[3] John P. Boyd, 1989: Chebyshev & Fourier spectral methods, Springer, 793 

pp. 



6 



Fig. 4. As in Fig. 3 but for t = 1.6037/5tt. 
[4] M.O. Deville, P.F. Fischer & E.H. Mund, 2002: High-order methods for 

incompressible fluid flow, Cambridge, 499 pp. 
[5] A. Fournier and D. Rosenberg, 2005: "Multiresolution adaptive spectral- 
element modeling of geophysical fluid dynamics using GASpAR," in prepa- 
ration. 

[6] Fournier, A., G. Beylkin and V. Cheruvu, 2005: Multiresolution adaptive 
space refinement in geophysical fluid dynamics simulation. Lecture Notes 
Comp. Sci. Eng., 41, 161-170. 

[7] Rosenberg, D., A. Fournier, P.F. Fischer & A. Pouquet, 2005: Geophysical- 
astrophysical spectral-element adaptive refinement (GASpAR): Object- 
oriented /i-adaptive code for geophysical fluid dynamics simulation. J. 
Comp. Phys., submitted. 



